Trace Formulas and Bogomolny's Transfer Operator 
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The trace formulas commonly used in discussing spectral properties of quantum or wave systems are 
derived simply and directly from the Bogomolny transfer operator. Special cases are the Gutzwiller 
formula, the Berry- Tabor formula, and the perturbed Berry- Tabor formula. A comment is made 
about interpolation formulas proposed in the latter case. 
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PACS number(s): 05.45.-a,03.65.Sq 
I. INTRODUCTION 

The Gutzwiller trace formulate and its relatives are the 
cornerstone of much of the subject of 'quantum chaology'. 
This formula gives a formal quasiclassical expression for 
the density of states, in other words, for the energy lev- 
els, of a,lhard' chaotic system. The Berry- Tabor trace 
formulaBu is the corresponding formula for an integrabLs 
system. These are descendants of the exact Selbergcl 
trace formula, which applies to spaces of negative cur- 
vature. 

Twenty-five years of work on the trace formulas leads 
to the conclusion that they are not especially useful for 
finding actual energy levels. This is because they are di- 
vergent formulas unless they are smoothed over energy. 
On the other hand, many interesting observables can be 
expressed in terms of the density of states smoothed over 
some energy scale. For example, the specific heat or 
magnetic susceptibility are observables of this type. The 
smoothing can come about because of finite temperature 
or finite experimental resolution, say. With appropriate 
smoothing, the trace formulas are very useful. 

We here provide a derivation of these and related for- 
mulas based on Bogomolny's transfer operatorQ. One of 
the main virtues of the transfer operator is that it can 
be used to study energy levels, wavefunctions and scat- 
tering amplitudes. The possibility of obtaining the trace 
formulas from the transfer operator is widely known, and 
certainly was to Bogomolny, especially in the Gutzwiller 
case. 

An explicit derivation of the Berry- Tabor formula or of 
the perturbed Berry- Tabor formula that is based on the 
Bogomolny operator does not seem to have appeared in 
print, however. That is the main purpose of this paper. 
We restrict consideration to the two-dimensional case, 
which is by far the most important in practice, and to a 
closed system with a discrete spectrum. _ 

The original derivation of the Gutzwiller formulaEI was 
based on a stationary phase evaluation of Feynman's path 
integral expression for the quantum propagator. This re- 
duces the propagator from a sum over all paths to a sum 
over just the classical paths. Although the Feynman for- 
mulation gives a persuasive picture of the relationship of 
quantum to classical mechanics, it does have some de- 
fects. For one thing, it is mathematically not very well 
defined. Also, its perspective on classical mechanics is 
very elementary. In particular, it does not incorporate 
the insights gained by use of Poincare's surface of section 
method. The present derivation removes these defects. 



The. .Berry- Tabor formula was derived in several 
waysBflfl. For example, the density of states expressed as 
a sum over explicitly known energy levels can be trans- 
formed by use of the Poisson sum formula. Action-angle 
variables are also useful to obtain this formulaB. 

With considerable generality, for a give* problem, that 
is, for a given Hamiltonian, there existsn a mathemati- 
cally well defined kernel or operator K(q,q',E), whose 
arguments q, q' are on a one-dimensional Poincare sur- 
face of section [SS] . The exact spectrum is determined 
by the Fredholm integral equation 



1>(q)= / dq'K{q,q',E)^(q'), 



(1) 



ss 



which has nontrivial solutions only if E is on the spec- 
trum. The integral runs over the surface of section. 
This is an exact quantum version of the classical SS 
method. A well known and much used examples is a 
billiard where the surface of section is the boundary and 
K(q, q' , E) is a Hankel function with argument kL{q 1 q'). 
Here k = \ / 2mE/h with m the mass and E the energy 
of the particle in the billiard, and L is the chord distance 
between points on the boundary specified by the g's. The 
exact or numerically approximated naethod is in this case 
called the boundary integral methooa. 



II. BOGOMOLNY'S TRANSFER OPERATOR, T 

In more general cases, K does not have a simple ex- 
pression. However, BogomolnjnJ in effect argued that in 
quasiclassical approximation K may be replaced by its 
asymptotic form for large separation of points q, q' . In 
this approximation, the operator, called T, is quite sim- 
ple, because it involves only a few, often just one, short 
orbits. Namely, 



T(q,q',E) = 



( 1 



\1~Kih 
x exp 



d 2 S(q,q',E) 



dqdq 1 



-S(q,q',E) + —v 



(2) 



where S is the action J pdq along the classical path of 
energy E, which leaves the surface of section at q' and 
returns to it for the first time at point q. 

Care must be taken with the phase attributed to the 
prefactor and with phase changes encountered at points 



along the path where the leading order semiclassical ap- 
proximation fails, e.g. at a billiard boundary. These are 
incorporated into the 'Maslov' phase v. It can be regarded 
as part of the action S, i.e. S — > S + hv/4 and thus is 
a quantum correction. Also, in doing integrals such as 
those below, contributions to v may be obtained. We 
shall, however, suppress v in the sequel, since that is not 
the point of this work. 

It should be noted that S is the generator of the clas- 
sical surface of section map, namely p = dS(q,q')/dq, 
p' = —dS/dq 1 , where p, p' are the momenta of the orbit 
parallel to the surface of section at q, q'. S thus con- 
tains all the long-time and long-orbit information of the 
classical mechanics. 



III. THE GENERAL QUASICLASSICAL TRACE 
FORMULA 



Eq. (Fy), with K replaced by T, has a solution when 
the Fredholm determinant!] 



D(E) = dct [1 - T(E)] = 0. 



(3) 



The denary of states is given by the logarithmic 
derivative 



d(E)- d(E) = d osc (E) = — Im 

7T 



d\nD(E + ie) 



dE 



> (4) 



where d(E) = J2 a $(E — E a ) an d d(E) is the smoothed 
(Weyl) density of states. Using the relationship lndet(l — 
T) = Trln(l — T), and expanding the logarithm as if T 
were small we have 



, ,„. 1 ^ 1 dr n (E) 



^ n dE 

n— 1 



(5) 



evaluated in stationary phase is given b 



>B 



T 2 (q,q')=J2 



1 



2nih 
i 



d 2 S 2p (q,q\E) 



x exp 



dqdq' 
S 2p (q,q',E) 



where 



S 2p (q, q', E) = S(q, q p ) + S(q p , q) 



(9) 



(10) 



is the action of a classical orbit from q' to q at energy E, 
which arrives at the surface of section at q for the second 
time after having crossed the first time at q p . The sta- 
tionary phase condition is d[S(q, q p ) + S(q p , q')]/dq p = 0. 
The solution of this equation, q p , is a function of q, q'. 
There are, in general, more than one solution q p to this 
stationary phase condition, which can lead to an expo- 
nential proliferation of orbits as a function of n. Thus, 
a T operator corresponding to multiple SS crossings can 
be defined, and it has exactly the same form as Eq. (0), 
although there are in general many terms of this type. 

There are two common situations. In the first, all n 
integrals may be performed by stationary phase. This is 
the Gutzwiller assumption. In the second case, the first 
n—1 integrals are well approximated by stationary phase 
but this method fails for the last integral. 

Eq. (0) may be iterated to obtain the final integral 



"n = / dq ^ 



1 



2Trih 



d 2 S p {q,q\E) 



l=l' 



x exp 



-S p (q,q,E) 



(11) 



where 



Here p denotes a 'closed' orbit from q, which crosses the 
SS n — 1 times and arrives back at q at its nth encounter 
with the SS. 



T n (E)=TvT(E) n . 



(6) 



Since T is in a quasiclassical sense close to a unitary 
operator, some of the eigenvalues of T are close to the 
unit circle, which makes the expansion (Eh divergent for 
real E. However, in this paper we ignore this problem, 
since, as we said, in practice the trace formulas are useful 
only when averages are taken, which suppress the contri- 
bution of the terms with large n. 

The trace formulas are obtained by evaluating the 
traces defining the r„'s in stationary phase approxima- 
tion, if this is meaningful. There are n integrals to be 
carried out, i.e. 



dqi . . . dq„T(q 1 ,q 2 ) . . . T(q n , qi). (7) 
The fundamental composition relation 

T 2 (q,q') = JdqT(q,q)T(q,q') (8) 



IV. THE GUTZWILLER TRACE FORMULA 

In the Gutzwiller case, the final q integral may be done 
by stationary phase also, that is, the phase S p (q, q, E)/% 
is a rapidly varying function of q. This picks out as 
stationary points the solutions of dS p (q,q' , E)/dq + 
dS p (q,q',E)/dq' = p - p' = at q = q' = q*. Thus 
periodic orbits are selected. Let q, q' be in the vicinity of 
q* and expand 

S p (q,q')~S p (E)+p*(6q-5q') 

+ \w 11 Sq 2 + W 12 SqSq' + ^W 22 5q' 2 , (12) 

where Sq = q—q*, 5q' = q'—q* andp* is the momentum of 
the periodic orbit. The W matrix depends on energy and 
on the orbit. Let the orbit p of 'length' n (i.e. returning 
to the surface of section n times) consist of r repetitions 
of 'length' s, with n — rs. The q integral is then done to 
obtain the contribution of this periodic orbit to the trace 



IF 



12 



Wn + 2W X2 + W 22 



exp 



- n S p{ E) 



(13) 



The factor s arises because the stationary point for the 
last integral can occur (or, in other words, the start- 
ing/ending point for the periodic orbit can be chosen) 
at any of s surface of section crossings. 

The prefactor is usually expressed in terms of the mon- 
odromy matrix M p of orbit p, which connects the final 
momentum-position p — p* — Sp — WnSq + W^Sq' and 
Sq to the initial Sp' = —W^Sq— \V22Sq' and Sq' . That is 



M„ 




Wf 2 -WnW22 
w\ 2 

w 12 



(14) 



so that the prefactor in Eq. 



1 



-1/2 



3h can be written 



|det(M p 

In the case of a repeated orbit, S p — rS s , and M p = 
Mg. In taking the derivative of t p only the rapidly vary- 
ing phase needs to be differentiated, giving a factor of 
dS s (E)/dE — T s , the period of the primitive orbit, lead- 
ing to the final expression 



Uosc — / 



Ts/fi 



|det(M;-l)| 



1/2 



cos 



S S (E) 



TV 

2 s 



(15) 



in which we restore the Maslov index. 



V. THE BERRY- TABOR FORMULA 

The above formula does not apply to integrable sys- 
tems. In that case, the periodic orbits are not isolated. 
The eigenvalues A of the monodromy matrix are, for iso- 
lated orbits, A = e ±7 , where for an unstable orbit the 
Lyapunov exponent 7 > 0, and for a stable orbit 7 is 
pure imaginary For an integrable system, 7 vanishes 
and formula (|l5|) is infinite. 

In the case of an intestable system, it is easiest to go 
to action-angle variablesa, 0,I,Q,J using as surface of 
section 9 = 0, and energy E = H(I, J) fixed. Then the 
action of the T operator becomes S(q,q') — * S(9 — 9'). 
The replacement of the action depending on two variables 
separately by a function of the difference of coordinates 
only characterizes the integrable system. In fact, the 
action is 



S{9-9',E) = {6-9')I + 2wJ. 



(16) 



The actions /, J specify the invariant torus, on which 
the orbit remains. The frequencies of rotation around 
this torus are 101 = dH/dl and ujj = dH/dJ . We 
may equally well use energy E and winding number 
u; = uji/ujj = (9 — 9')/2tt as variables to specify the orbit. 
[This fails for harmonic oscillators where the frequencies 
are constants. Modifications are necessary in that caseo.] 
Thus, I, J are regarded as functions of E, 9 — 9' . It is 
easy to see that dS/d9 = /. Obviously, the eigenfunction 
of the T operator is in this case ip — expiI9. 



The first n — 1 integrals may be carried out by sta- 
tionary phase, as before, giving a result S p (6 — 9') where 
6, which is physically equivalent to 9', is set equal to 
9' + 27T771. The winding number is m/n, i.e. it is rational. 
Clearly S p (9 - 9') = nS(A9), where A9 = (9 - 9')/n. 
Remark first that in this limit p is a periodic orbit. It 
is a periodic orbit for each 9 with identical period and 
other properties, in other words, this describes a set of 
periodic orbits, which are not isolated. 

The final integral Eq. ( pT| ) is then trivial, giving a fac- 
tor 2ir. Note a factor i _ 5 , which remains and gives a shift 
7r/4 in the final answer below. The prefactor involves the 
second derivative (at constant E), S" — -^dI/dA9. This 
has traditionally been expressed in terms of the function 
9e(I) = J, which determines the action J given the en- 
ergy and /. One finds that g' E = —A9/2tt, the derivative 
being taken at fixed E. Therefore S!l = —(2img'g)z} n 
Putting this all together we obtain the standard resultn-El 



&osc — / 



T 



Trh 3 ' 2 n 3 / 2 \g%\ 1/2 



cos 



o n 



IT 7T 

2 p 4 



(17) 



There are other cases, in which the last integral Eq. 
(O) cannot be done by stationary phase. There may 
be particular sets of nonisolated orbits in an otherwise 
chaotic system. Prominent examples are the bouncing 
ball orbits in the stadiumllj or Sinai billiardEi These are 
also nonisolated orbits in the ray splitting billiardlij. 



VI. PERTURBED BERRY- TABOR FORMULA 

Another important example is the perturbed inte- 
grable system. Then the action defining the T operator 
has the form S(9 - 9') + eWi{6,0'), and e « 1. To the 
leading order in e, we may assume that the stationary 
phase points of the first n — 1 integrals are not shifted. 
The last integral then involves 



I w = — I d9 exp 



£* W 



(18) 



where W(9) = Er=i^i(^ - r -i,6 - r ), and 9 r = 
2irmr/n for the winding number m/n. It can be shown 
that W{9) = W„{9 1 9 - 2ttto), where W n {9,9') is the 
first order correction to the action for the orbit that re- 
turns to the SS for the nth time. The perturbed re- 
sult is just the Berry- Tabor result with the substitution 
cos <p — ► Re [Iyy ps(;p(*0)]j where <f> is the argument of co- 
sine in Eq. @E3. 

Clearly, Iw is a sort of generalized Bessel function. 
Formulas for Iw interpolating between e = and e/h 
large have been given in the literature. The first were 
based on the idea that often W will be_Hiell approximated 
by the first terms of its Fourier seriesEj, e.g. W ~ wq + 
wiCOs(9 — 0o)- The phase shift 9q is of no importance, 
and Iw ~ ex.p(iewo/h)Jo(ewi/h), where Jo is the Bessel 
function. 



Interpolation formulas along this line using more pa- 
rameters have been given in the literatureli-l. We give 
an alternative four parameter formula in the Appendix, 
since we believe the existing formulas are defective. 

One comment concerning these formulas based on a 
small parameter e is in order. Both classical perturba- 
tion theory and quantum results based on small e have 
relatively large effects, proportional to y/e, when the sys- 
tem is in the neighborhood of a resonance. Furthermore, 
the trace formulas apply precisely to the resonant cases. 

Nevertheless, it is the case that for the perturbed trace 
formula, the corrections incorporated into Eq. (flq) can 
be expanded in powers of e, and such a series is conver- 
gent. 

The reason is as follows: To find energy levels and 
wave functions, one must solve Eqs. (H,p]). This, in ef- 
fect, depends on the very high terms and convergence 
properties of the trace formulas. The. -y/e appears if one 
is interested in this level of accuracycJ. However, to look 
at long range energy level correlations, and the smoothed 
density of states, it is only necessary to consider the lead- 
ing terms in the trace formula expansion. These leading 
terms can be calculated from short periodic orbits. Al- 
ready classically, it is known that short time quantities 
can be calculated perturbatively as a power series in e. 



second derivative, evaluated at an extremum, is W" — 
-(£') 2 WW cos £. Since d0/(%= l-Acos£-2.Bcos2£, we 
find C(0) = [1 - A - 2B]~ l and £'(tt) = [1 + A - 2B]' 1 . 
We may solve for A, B, W^°\ and W^ to obtain a 
parametrized version of W , which interpolates between 
e/% = and e/% large. 

The integral Eq. (jl^) now is 

Iw = -^- [ d£(l- Acos£-2Bcos2£) 

2tt J 

xexp[i(e/?i)VF (1) cosd , (Al) 

where Iq = exp \i (e/%) W 1 - ^ . The integrals may be ex- 
pressed as Bessel functions J m , whose arguments are 
eW^'/h, giving 



iw 



I (J -iAJ 1 +2BJ 2 ). 



(A2) 



As an exercise, and to compare with UGT, we calcu- 
late Iw for coupled quartic oscillators. The Hamiltonian 
for this case is 



H 



1 



(pl + pl) + ax A + by 4 + ex 2 y 2 



(A3) 



VII. CONCLUSION 

The Bogomolny transfer operator T provides a unified 
approach to the quasiclassical description of wave sys- 
tems, without needing to make specific reference to the 
chaotic or intcgrable properties of the system. Its great- 
est importance is, no doubt, that it provides a method to 
find actual energy levels and wave functions. However, 
it can also be used to find the trace formulas, which are 
widely used in this field. This derivation of the trace for- 
mulas is quite easy and direct and it is on at least as firm 
a mathematical footing as most other derivations. 



APPENDIX A: UNIFORM APPROXIMATION 

Suppose W(9) has a single maximum and minimum 
at 9 = 9 x ,9 n , respectively. For large e/% the integral 
( |l8| ) can be done by stationary phase. This is sometimes 
called the Gutzwiller case because the continuum of iden- 
tical periodic orbits has turned into two isolated orbits. 
The integral in this approximation depends on four pa- 
rameters, namely, the magnitudes of W and W" at the 
extrema. It does not depend directly on the position of 
the extrema. Thus, in order to have the correct large e/% 
limit, one should use-a four parameter expression for Iw- 

Inspired by Ref.li-3, [UGT], we parametrize W(9) = 
W(0) + W (i) cos[£(0)], where £($ x ) = and f(0 n ) = n. 
Note that W^ and W^ are not necessarily equal to 
Wo, w\- We may as well take 9 X = and 9 n = n, 
to simplify things. Now take £ as the solution of 6 = 
£ — Asin£ — i?sin2£. This goes beyond UGT, who take 
B = 0. The derivative is W' = -£'(6>)W (1) sin£, and the 



The correspondence with UGT goes as follows: a ^ 
a(X)/b, b ^ a(X)b, e ^ 2Aa(A), where a(A) is a speci- 
fied function of A. However, we may choose the unit of 
length such that a(X) = 1, i.e. a = 1/6. Thus A, the small 
parameter of UGP, is in effect our e/2. 

Using action-angle variables, with A9 = 2ir, we find 



S(6 - 9') = 



4KE^ 4 
3^r 



(2n) 4 



-,1/4 



(A4) 



where K = K(rri-^= 1/2) is the complete elliptic integral 
of the first kindt£l. Here E is the energy of the unper- 
turbed system. The perturbing action W n is found (to 
leading order in e) in terms of Jacobi elliptic functionsE£l 
as 



W n (9,9 , ) = 



u> 



KE 3 / 4 



(A5) 



-1-KThUJ 



d(9"sd 5 



[a{9" + 9)}sd 2 [a9"/uj] 



where a = 2K/ir and to — (9 — 9') /27m. We focus on the 
1-1 resonance, w = 1, and W{9) = Wi(0,9 - 2n). Note 
that the shape of W(9) for fixed ui is independent of a 
and E and depends only on the winding number. 

Fig. p] gives a graph of W(9). It has two maxima and 
minima, which are symmetrically related. It is easy to 
extend the interpolation expression to take this into ac- 
count. Re[/iy] as a function of e/% is presented in Fig. g, 
and Im[iV] is shown in Fig. pi We find that the inter- 
polation formula and the exact integral agree to within 
a part in 1CU 5 (inset Fig. 0). 




FIG. 1. The function W(9) vs. 9 for anharmonic oscilla- 
tors with x 2 y 2 coupling, for the (1,1) resonance. The param- 
eters a, b, and E are defined in the text and each is taken to 
be unity. 

UGT assume B is zero. This means that they cannot 
generally satisfy the asymptotic condition for large e/h. 
It can be seen that B = implies a relation between the 
curvatures at minimum and maximum, namely, 



vV«/^' lin +yW) 



W" 

r r max 



= 2/r, 



(A6) 



where r is the number of maxima (or minima) of W(9) 
per 2-7T (r = 2 for 1-1 resonance). However, for quartic 
oscillators with x 2 y 2 coupling, this is well satisfied, the 
l.h.s. of Eq. (|A6|) being equal to 1.00021. This probably 
accounts for the good numerical results of UGT in this 
case. We find A = -1.49 x 10~ 2 , B = -1.04 x 10~ 4 , 
independent of the other parameters, and W^ = —0.46 
and WW = 0.27 for E = a = 1/6 = 1. 



depending on the invariant tori, defined only at e = 0. 
Its definition cannot be easily extended for e > 0, since 
the topology of the invariant tori changes under pertur- 
bation. We found the simple formula, in the present ex- 
ercise, 



9e{I) = J 



fr-©"> 3 



3/4 



(A7) 



UGT do not mention the curvatures W":„, W[' 

linn 7 iiidx 



explic- 
itly. These parameters are replaced by the determinants 
of the monodromy matrix Mp — 1, at the stable and un- 
stable orbits [see our Eq. dig )]. These determinants can 
be expressed in terms of W" and g E . Eliminating thus 
W", we obtain the result of UGT in our notation, 



det(M„ - 1) = 
det(M s - 1) = 



-ceg'E 



ceg's 



W {1) /[l-A + 2Bf 
W {1) /[l + A + 2B] 2 . 



(A8) 



Here c is a known constant. 
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FIG. 3. The imaginary part of Iw as a function of e/h. 

Given g'^ and W^- 1 ' , we can find A and B in terms of 
the determinants. Setting B = leaves too few parame- 
ters, so g'g, extended to be a function of e, is pressed into 
service by UGT. However, for small e, the determinants 
are proportional to e because the Lyapunov exponents 
are proportional to y/e. Their ratio is a constant not nec- 
essarily equal to —1. The explicit e cancels out leaving 
g^(e) as nearly constant for small e, but not necessarily 
close to the zero e limit, unless the system is such that 
B happens to be small. 



10 
e/h 



FIG. 2. The real part of Iw as a function of e/h. The differ- 
ence between this function computed by numerical quadrature 
(Iw) an d from the interpolation formula (I w ) is shown in the 
inset. 

To obtain a fourth parameter, UGT formally take g'^ 
as an empirical parameter depending on e, arguing that 
an independent evaluation of gE is rather laborious and 
time consuming. However, g^ is a geometrical quantity 
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